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Abstract 

We solve the first-passage problem for the Heston random diffusion model. We obtain exact 
analytical expressions for the survival and hitting probabilities to a given level of return. We 
study several asymptotic behaviors and obtain approximate forms of these probabilities which 
prove, among other interesting properties, the non-existence of a mean first-passage time. One 
significant result is the evidence of extreme deviations -which implies a high risk of default- when 
certain dimensionless parameter, related to the strength of the volatility fiuctuations, increases. 
We believe that this may provide an effective tool for risk control which can be readily applicable 
to real markets. 
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I. INTRODUCTION 



The study of first-passage and exit problems liave a long and standing tradition in physics, 



engineering and natural sciences [l, Perhaps the most important example of an exit 
problem in physics is provided by the "Kramers problem" where one studies the possible 
escape, owing to noise, of a system from a stable estate 4]. Classical examples of first- 
passage problems are the collapse of mechanical structures because of random external 
vibrations which attain an extreme amplitude beyond the stability threshold; or the "false 
alarm problem" where internal fluctuations induce the current or voltage of an electric circuit 
to reach a critical value for which an alarm is triggered j^, [g]. 

In finance, the study of extreme events from the perspective of first-passage or exit 
problems is still beginning. In few recent papers we have address the escape problem in the 
context of econophysics using several techniques for different settings B,l3y,y,Q. In the 



present paper we want to address for financial series the related issue of the first-passage 
problem. This is a relevant question because it is intimately connected to extreme values 
and, therefore, to risk and default. 

Needless to say that the evaluation and subsequent risk control should be central issues 
in finance. Crisis are, nevertheless, inherent characteristics of financial markets although 
one expects them to be rare. Unfortunately this trend is changing, at least in the last two 
decades, for the frequency of appearance of extreme events seems to be increasing and there 



is a growing consensus that there shou 
option pricing and structured products 



d be a change of model, specially in the fields of 



12|. 



Traditionally tasks related to risk control have been based on the assumption that prices 
follow the geometric Brownian motion (GBM). However, it is known for long 13|, that 
the Gaussian distributions resulting from the GBM decay for large values of the price (or the 
return) much more steeply than the empirical distributions. In other words, real financial 
distributions present fat tails and prices are highly leptokurtic 15| . 

Therefore, the GBM model clearly underestimates risk. Indeed, risk is associated with 
large deviations of returns and the appearance of these deviations is determined by the 
tails of the probability distribution. Since for the GBM tails are much thinner than the 
empirical ones, we conclude that actual risk is higher than the risk foreseen by the GBM. 
Practitioners usually remedy this (dangerous [12;]) inadequacy of the GBM model by the 
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ad-hoc procedure of assuming that prices follow a mixture of two Gaussian distributions. 
The mixture is tailored to reproduce the real probability distribution, that is to say, one 
Gaussian curve adjusts the center of the real distribution and the second one the tails 16|. 
Whereas this certainly improves risk management, it is an unsatisfactory procedure, for it 
is little more than curve fitting without any model sustaining it. Similar objections can be 
raised to other approaches based on adjusting actual prices by truncated Levy distributions 



171, 1181, iJi 

Another facet of the problem is provided by the fact that the GBM assumes a constant 
volatility. The latter defined as the standard deviation of returns. However, the recording 
of the empirical prices of financial options clearly indicates that the volatility im plie d in 
these prices is not a constant, not even a function of time, but a random variable 20|, l21 |. 
This gave rise in the late ninety eighties to the so-called stochastic volatility (SV) models 
for which return and volatility constitute a two-dimensional diffusion process governed by a 



pair of Langevin equations 22, l23|, |24 1 . 



Financial models based on SV have been gaining increasing acceptance, not only for 
the initial motivation of correcting option prices, but because they provide a natural and 
reasonable explanation to many empirical observations which are gathered together under 
the collective name of "stylized facts" In particular SV models result in fat-tailed 

distributions of returns which overcomes one of the main objections against the GBM. 
We should also note that volatility is today a crucial concept in any financial setting and 
many financial products are based on it. Therefore, the risk associated with volatility is 
particularly important from theoretical, as well as practical, aspects of the problem. 

The measure of risk is provided by the probability of attaining an extreme value (or a 
preassigned critical label) and, in the standard approach, this probability is obtained through 
the distribution of prices. In such a procedure the probability of reaching certain critical 
mark L is secured by counting all times the process has reached the label L. There is an 
alternative (and finer) approach to the problem which consists in counting only the event 
of reaching L for the first time. This is the first-passage problem that, as we have explained 
above, is an old and well known acquaintance of physics, engineering and natural sciences 

m 

. Within the context of continuous-time random walks, we have recently shown that 
risk control methods based on first-passage techniques are more efficient than traditional 
ones, since the latter clearly underestimate risk 10|. 
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Herein we address the first-passage (or hitting) problem for financial time series under 
the assumption of stochastic volatility. We choose one particular SV model, the Heston 
model, because it has the advantage of allowing exact analytical results 
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. In 



a very recent paper [11] we have solved a closely related issue: the escape problem for the 
Heston model. In that problem the process can exit, for the first time, a given interval and 
this leads to a two-barrier problem. However, in the first-passage problem the process only 
crosses a given critical value which amounts to solving an one-barrier problem. 

To our knowledge there are a few but increasing number of works dealing with first- 



passage problems [29|, |30|, |3ll] and exit problems 



10 



11 



32 



33 



34| in finance. Many 



of them address the exit problem in a restricted way, because they only obtain the mean 
exit time, i.e., the first moment of the exit distribution. This moment is easier to handle 
than the entire distribution, although it provides much less information, specially about the 
time evolution of the exit problem. On the other hand first-passage problems are usually 
more intricate because in many cases the mean first-passage time does not exist. This is the 
case, for instance, of the Brownian motion or the standard SV models (see below). In such 
situations one is either bound to get the entire first-passage distribution, which is usually 
quite involved, or to introduce ad-hoc forces for securing the existence of the first moment 
of the hitting 
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33|. 



In this paper we solve the complete first-passage problem for the Heston SV model. We 
obtain the exact expression of the survival probability and show that the mean first-passage 
time does not exist. We also average the volatility in order to work out the problem for 
the return alone. We obtain approximate expressions of the survival probability in several 
asymptotic regimes -long times, large volatilities and small volatility fluctuations- and find 
that in these cases the (approximate) survival probability is Gaussian and has the same 
form than that of the Wiener process. We also obtain the asymptotic expression for survival 
when the fluctuations of the volatility are high. In the later case the Gaussian character is 
lost and the probability of hitting an extreme value is much higher than in the case of mild 
volatility fluctuations. 

The paper in organized as follows. In Sec. [Ill we present the Heston model and obtain 
the exact solution to the first-passage problem. In Sec. IIIII we obtain the approximate 
expression for the survival probability in several asymptotic limits. In Sec. IIVI we treat 
the problem according to the fluctuations of the volatility. In Sec. |V] we average out the 
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volatility assuming it has reached the stationary state. This allows us to get exact as well as 
approximate expressions for the survival probability of the return alone. A brief summary 
of the main results along with few conclusions are given in Sec. I VII Some more technical 
details are in an Appendix. 



II. THE HESTON MODEL AND THE FIRST-PASSAGE TIME DISTRIBUTION 



Let X{t) be the zero- mean return X{t) defined through the stochastic differential (in the 
Ito sense): 

-«^^-{^>' « 

where P{t) is a speculative price or the value of an index and (■) denotes the average. In 
terms of X{t) the Heston model is a two-dimensional diffusion process {X{t),Y{t)) 
described by the following pair of stochastic differential equations (again, in the Ito sense) 

dX{t) = ^Y{t)dWi{t), (2) 



dY{t) = -a [Y{t) - m'\ dt + k,jY{t)dW2{t), (3) 

where Wi(t) are Wiener processes, i.e. dWi{t) = ^i(t)dt(i = 1,2), where ^i(t) are zero-mean 
Gaussian white noises with {^i{t)^i{t')) = 6ij6{t — t') |35 |. 

Equation ([2]) shows that Y{t) is the variance of return and the volatility a(t) is given by 



a{t) = ^Y{t). (4) 

However, as long as no confusion arises, we will use the term "volatility variable" or just 
"volatility" for the random process Y{t). Moreover, as proved many years ago 36|, the 
volatility process Y(t) defined by Eq. is positive and the volatility is well defined. 

In Eq. ([3]) the parameter m is the so-called "normal level of volatility", a > is related 
to the "reverting force" toward the normal level m (which cannot be zero, see discussion 
below) and k, sometimes referred to as the "vol of vol", measures the fluctuations of the 
volatility. It is useful to keep in mind that the quantities Y{t), a, m^, and k have them all 
units of 1/times; while the zero-mean return and, consequently L are dimensionless. 

One can easily see that Y(t) is a homogeneous and stationary random process whose 
probability density function as t — > oo (or, equivalently as to —oo, where is the initial 
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time) is given by the Gamma distribution: 

pM = ^^^y^-e-(-/'^^)^ (5) 

where u = (2a/A;^)m^. 

From Eq. ([5]) we readily see that the stationary mean value and the stationary variance 
of the volatility variable are 

(F(t))st = VaT{Y{t)}^^ = I2a)m\ (6) 

As long as Y{t) is a (positive) random variable neither its variance nor its mean can be zero, 
hence m 7^ 0. 

After this brief summary on the Heston model we turn our attention to first-passage 
problems. Let S{x,y,t) be the probability that the zero mean return, X{t), initially at 
X{0) = X js^l with volatility 1^(0) = y, has never crossed the critical level L before time t. 

When L < X (i.e., L is a level of losses) S{x,y,t) coincides with the survival probability 
(SP) for the joint process {X{t),Y(t)) to be, at time t, still inside the semi-infinite strip 

L < X{t) < 00, < Y{t) < 00. 

In such a case the (first) hitting probability, that is, the probability of having a loss labeled 
by L (if L ^ —00, this is the "default probability") is given by 

W{x,y,t) = l-S{x,y,t). (7) 

When L > X, S{x, y, t) is the SP of the semi-infinite strip 

-00 < X{t) < L, < Y{t) < 00, 

and the probability of hitting a profit L is also given by Eq. (J?!). Since in our Heston model 
we do not consider any bias both situations -loss and profit- are symmetrical. 



The survival 
value problem 



probability S{x, y, t) is the solution to the following initial and boundary 



dS , on-^S 1,2 d'^S 1 



with initial and boundary conditions respectively given by 



S{x,y,0) = l, S{L,y,t) = 0. 
6 



(9) 



Note that the symmetry between losses and profits just mentioned can be readily seen by 
the fact that under the change of variable 



z=\L-x\ (10) 



problem (IHI)-© remains unchanged: 



with 

S{z,y,0) = l, and S{0,y,t) = 0. (12) 
Since z >0 and S{0,y,t) = 0, we can define the Fourier transform 



oo 



S{uj,y,t) = / S{z,y,t) sinujzdz, (13) 





which turns Eqs. (|TT]) - (|T2l) into the following initial value problem 
with initial condition 

S{u,y,0) = VilM, (15) 

where 

POO 

V{l/uj)=j^ sinujzdz = l/uj (^7^0), (16) 

is the Cauchy principal value js^. 

After the definition of the dimensionless variables 

T = at, V = y/a, (17) 

the problem above reads 

|^ = _(„_9,^ + (^V2)0-(-V2H.S. (18) 

and 

S{u,v,0)=Vil/cu), (19) 

where 

e = m'^/a, (20) 
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is the (dimensionless) normal level of volatility and 



P = k/a, 



(21) 



is a dimensionless parameter which gauges the volatility fluctuations with respect to the 
deterministic strength, measured by a, toward the normal level. 

As can be easily seen by direct substitution the solution to the problem posed by Eqs. 
(HHD-din]) is furnished by 



where 



S{uj, v,t)=V (I/lu) exp {-A{u, r) - 2B{u, t)v/P^} 



A{iu,T) = {29/ /3^) / B{uo,s)ds, 
Jo 



and B{uj,t) obeys the Riccati equation 



B = -B-B^ + {j3u/2f 



(22) 



(23) 



(24) 



with initial condition i?(co',0) = 0. 

The explicit expression for the function A{uj, r) and -B(a;, r) are obtained in Appendix 
lAl Thus 

"/i+(cj) +/i_(cu)e^^(-)"" 



A{uj,t) = (20/^) /i_(cu)r + ln 



Muj) 



and 



where 



I _ e-A(w)7 



Hence (cf Eqs. (l22l). (|25|) and (|26l)) 



l + [/i_(^)//i+M]e-^(-)-' 

/i±(a;) = [A(^)±l]/2. 



(25) 



(26) 



(27) 



A(cu)e-'^-('^)^ 



26I//32 



(28) 



/i+(u;) + yU„(a;)e~^H^ 

The solution to the hitting problem for the two-dimensional Heston SV model is therefore 
given by the Fourier inversion 

2 ~ 

S{z,v,t) = — S {u, V, t) sin uzdu. 
IT Jo 
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FIG. 1: (Color online) Log-log representation of the survival probability S{z,v,t) as a function 
of z for T = 0.5 (around 11 trading days), v = 9 and /? = 0.1. Empty diamonds represent the 
numerical evaluation of the exact expression gien in Eq. (|30p . Solid circles represent simulation 



2g] which 



points of the Heston model. We use the realistic parameters values presented in Ref. 
are similar to those encountered in actual markets: a = 0.045 1/day, = 8.62 x 10~^ 1/day and 
hence 6 = 1.92 x 10"^. 



Plugging into this equation the expression for S{uj,v,t) given in Eq. (!28|) and taking into 
account Eq. ( fT6l) . we get 



S{z,v,t) = - f 

TT JO 



°° smuz 



UJ 



(29) 



ix+{uj) + //_(a;)e-'^('^)^ 

In Fig. [1] we represent the numerical evaluation of this Fourier integral. We also show 
there the SP obtained through numerical simulations of the Heston model 40|. The match 
between both results is an indication of the correctness of our calculations. 



III. ASYMPTOTIC APPROXIMATIONS 



Equation (l29l) provides the complete solution to the first-passage problem for the Heston 
model which can be used to get numerical values for the SP and hitting probabilities as is 
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shown in Fig. [T] as well as in the rest of figures of the paper. However, it is not the most 
convenient form of S{z, v, r) for bringing into light general properties of the first-passage 
problem. We will, therefore, try to find approximate expressions that may uncover as many 
features of the SP as possible. We first obtain asymptotic expressions of the SP for long 
times which, in turn, will prove the non existence of an average hitting (i.e., first-passage) 
time. Subsequently we will obtain approximate expressions of the SP valid for large to 
moderate volatilities. 

We should also note that both approximations -i.e., long times and large volatilities- 
coincide. As we will see next this coincidence is due to the particular form of the functions 
/i_(ci;) and B{uj,t) defined above. 



A. Long-time asymptotic expressions and mean-first passage time 

We now focus on the behavior of the SP in the asymptotic limit when t — > oo, to this 
end rewrite Eq. (!29|) in the following form: 



^, , 2 r°° smuz 

S{Z,V,T) = - 



n Jo (jj 



For times such that 



A(tu)e-^('^'^)''/^ 



{9/P')t » 1, 



261//32 



-(2er//32)M-H^^_ (30) 



we are entitled to use the saddle-point method to secure an approximation of the integral 
in Eq. fl5U]) . One can easily show that = is a minimum of //-(cu) (cf. Eq. (1271) ) and 
expanding around this minimum yields 



fx4uj) = {Puj/2f + 0{uj^). 



Then 



S{z,v,t) 



2 /■°° sincu^; 



TT JO id 



-A{uj)t 



261//32 



(31) 



(32) 



Following the saddle-point method 4l| we expand to low orders in u the argument of the 
exponential term (see Eq. (!26|) ) 

B{u;, r) = {f3uj/2)\l - e"") + 0{u^). (33) 

Then taking into account Eq. fl3T]) and (cf. Eq. fl27|) ) 

A{iu) = l + {l3uy/2 + 0{uj^), /x+(cu) = l + (/5cu/2)2 + 0(a;^), (34) 
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we have 



AM 



-{2/f3'^)B{ui,v)v 



-1 29/13^ 



Substituting Eq. ([35]) into Eq. ([32]) yields 

2 sin uz 



1 + 



exp 



b[z,v^T) — — I exp 

TT JO UJ 



(35) 



(36) 



where 



\{T,v)=2eT + 2{l-e-^)v. (37) 

I — I 

The integral appearing in the right hand side of Eq. (1361) can be performed exactly [42] and 



we get 



where 



S{z, V, r) ~ Erf 



\iv,T) 
2 



is the error function 



43|. 



Erf (x) = — f 

x/TT Jo 



(38) 



It is interesting to note that this approximation to the survival probability fulfills not only 
the boundary condition but the initial condition as well. The latter being quite remarkable 
since the approximation is based on the assumption of long time. That Eq. (138]) obeys the 
boundary condition is readily seen from the fact that Erf(O) = 0. As to the initial condition 
we have S{z,v,t = 0) = Erf(oo) = 1. All of this is an indication of the goodness of the 
approximation as is corraborated by Fig. [21 

We shall now look at the mean first-passage time. Let us denote by ti^^p the time when 
the return X{t), starting from a known value x 37|, attains a certain (critical) level L 
for the first time. This time, usually referred to as the first-passage (or first-hitting) time, 
depends on x, L and the current level of volatility y. It is evidently a random variable 
whose average T = (tpp) is termed as the mean first-passage time (MFPT). In terms of the 
survival probability this average hitting time is given by 38 ] 



T{x,y) 



S{x, y, t)dt. 



Mathematical analysis tells us that in order for this integral to exist it is necessary that 
S{x,y,t) decay faster than 1/t as t ^ oo. Now looking at Eq. (l38l) and recalling that 
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Erf(x) = 2x/^ + 0{x^) we see that the SP of the Heston model decays when r ^ oo as 

S{z,v,t)^-^[1 + 0{1/t)]. 



Therefore, the survival probability falls off as 1/vt and the MFPT, 

T{x,y) = oo, 

does not exist. 

B. The first-passage problem for moderate and large volatilities 

Let us start by observing, as it appears in the right hand side of Eq. fl29l) . that the 
dependence of the SP on the volatility is that of a linear exponential. This suggests, as in the 
previous case when {9 / (3'^)t ^ 1, that as long as v is not small -i.e., from moderate to large 
volatilities- we may utilize again the saddle-point method for performing an approximate 
evaluation of Eq. (12^ . 

The case of large volatility turns out to be completely analogous to the long time approx- 
imation just discussed. In fact both cases lead to exactly the same approximate expression 
for the SP. Indeed, from Eq. fl33l) we see that, like fi^{uj), B{t,v) also attains a minimum 
value at = 0. Then, if v is large enough, the saddle-point method allows us to write the 
exact SP, Eq. (l29l) . in the approximate form 



, 2 r°°smujz 

S{Z,V,T, 



TC Jo LU 



^_,2(i_,-.)„/2^^_ (39) 



As before we expand to the lowest order in uj the bracketed expression in the integrand; we 
have (see Eqs. (1^ and (IMl)) 



^2, 

g-(/3c^/2)2r 



1 + O {u') 



H+{u) + /i„(c<;)e-^HH- 1 + {Puj/2ye-' 
Substituting this into Eq. (15^ we obtain again Eq. (15^ . which proves that the approximate 
expression of the SP valid for moderate to large volatilities is also given by Eq. (!38|) : 

^(^,t;,r) - Erf [ ^ ) , (t; > 1). (40) 

The correctness of this approximation is clearly shown in Fig. [31 
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FIG. 2: (Color online) Log-log representation of the survival probability as a function of the scaled 
time r with z = 0.01 and (3=1. Circles, squares and diamonds represent the numerical evaluation 
of the exact expression, Eq. (j30p . for the three values of the initial volatility v shown in the figure. 
Solid lines correspond to the asymptotic expression given by Eq. (|38p in terms of the error function 
and different curves correspond to different values of v. We see that as the volatility becomes 
smaller the exact expression needs a longer time to converge to the asymptotic expression. When 
the initial volatility is five times larger than the normal level, the convergence time is about 10^ 
trading days. For v/6 = 100 this time is much less (about 10^ days). Let us note that for volatilities 
thousand times higher than the normal level, the asymptotic curve follows very closely the exact 
SP for all values of r which is in agreement with the results of Sec. IIIIBl We use the realistic 
parameter values given in Fig. [TJ 

In the original units (cf. Eqs. (ITU]) and (IT7|) ) the approximate SP reads 



We remark that this approximation is valid for both long times and large volatilities. In the 




(41) 



where (cf. Eq. (EZ])) 




(42) 
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original units the condition of large volatilities, f ^ 1, is simply 



1/ > a; 



while the condition of long times, {9 / (3'^)t ^ 1, reads 

2 



t > 



k 



am 



that is, the time must be longer than a characteristic time, t^. = {k/amy, formed out of the 
three parameters of the Heston model: the vol of vol k, the reversion to the mean a, and 
the normal level m. 

As long as we assume that neither the volatility nor the time are small, we can safely 
use approximation fHOj) to obtain the behavior of the SP as \L\ —>■ oo. This is, perhaps, 
one of the most interesting aspect of the first-passage problem to look at, since the limit 
\L\ oo is intimately related to extreme risks and, ultimately, to default. Using the 
following asymptotic expression of the error function j43 | 



Erf(x) 



we have 



S{z,v,t) 



z0r 



l + 0{l/x^ 



(43) 



The hitting probability to the level L, W{z,v,t), is given by Eq. ([7]). The default or 
uprising probabilities will be obtained by assuming |L| — oo. In this way we find the 
following exponential decrease 



w 



i\L\ 



oo . 



(44) 



We finally take a closer look at the behavior of the SP with increasing values of the 
volatility. Equation ( l40l) is suitable for this purpose because it is valid for large values of 
the volatility. Intuition tells us that as ^ oo the survival probability should tend to zero 
which is readily seen from Eqs. fl37j) and fj40l) . Let us elucidate how this limit works. We 
first note that 



1 + 



4(1-6- 



+ Oil/v 
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FIG. 3: (Color online) Log-log representation of the survival probability S{z,v,t) as a function of 
the initial volatility with z = 0.01, and (3=1. Diamonds and squares correspond to the numerical 
evaluation of the exact expression (j30p at two instants of time r = 0.1 and t = 1, 000. Solid lines 
represent the asymptotic error function given by Eq. (j40p . Note that the approximation improves 
as time increases, in agreement with the results of Sec. IIII A[ We also observe a power law decay 
of the form with increasing volatility in agreement with the asymptotic estimate given by 

Eq. (|15|) . We have taken the realistic parameters of Fig. [TJ 



Moreover 43 1 



Erf(x) = ^[x- xV3 + 0{x% 



Merging these two expressions into Eq. f HOl) we get 



2z 

S{z, V, t) ^ , ^ [1 + (l/v)] . (45) 

'nil — e~'^)v 



Therefore, when volatility v increases the SP decreases as v ^1'^ . Figure [3] shows this power 
law behavior. 
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IV. THE FIRST-PASSAGE PROBLEM ACCORDING TO VOLATILITY FLUC- 
TUATIONS 

In the last section we have been deahng with approximate solutions to the hitting problem 
depending on the asymptotic values of time or volatility. We shall now look into other 
interesting and useful approximations. 

Recall that the parameter 

/3 = k/a 

measures the strength of the volatility fluctuations -given by the vol of vol k- in relation 
to the deterministic pull -measured by a- toward the normal level of volatility. In this way 
we may term the dimensionless parameter (3 as the "normalized volatility fluctuation" . On 
the other hand from Eqs. ([6]) and fl2T]) we see that (3'^ = (2/a;m^)Var{y(t)}st. In other 
words, is proportional to the amplitude of the volatility autocorrelation in the stationary 
state. Therefore, the higher (3 is, the more intense the stationary autocorrelation becomes. 
We, therefore, expect two different patterns according to whether (3 is small or large. This 
duality is also present in the stationary distribution of the volatility. In effect, we see from 
Eq. fl55l) below that the behavior of Pst{v) drastically changes as (3 goes from small values 
bounded by /5 < (26*)^/^ to larger values such that (3 > (26')^/^. Indeed, as v ^ the former 
case yields Pst{v) —>■ while in the latter Pst{v) oo; a fact that has to be taken into 
account in the numerical simulations of the Heston model [4^. 

We will explore which are the consequences of this on the first-passage problem of the 
return X(t). Let us first investigate the case /3 — > 0. 

A. Small fluctuations 

We start form the exact expression of the SP given in Eq. ( I29l) and suppose that volatility 
fluctuations are weak. Specifically, we assume that the vol of vol is much smaller than the 
pulling toward normal level, k <^ a; that is 

/9< 1. 

A glance at Eq. (1271) suffices to realize that the assumption of /3 — >■ is equivalent to 
assuming uj 0. But the latter is precisely the limit taken in the previous section as a 
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consequence of the saddle-point approximation. We therefore expect that the case of weak 
fluctuations will lead to the same approximation than those of long times or large volatilities. 
Let us briefly show that this is indeed the case. 

In effect, we start from Eq. fl2I?l) which we rewrite in the form 

-1 26»//32 



S{z,v,t) 



2 sin uz 



TT JO 



UJ 



X exp I 



[efx^{uj)T + B{uj,r)v]}duj. 



(46) 



If /? — the integral can be approximately performed by the saddle-point method which 
implies the expansion of the exponential term around its maximum located as before at 
a; = 0. In other words: 

exp{- [ei24uj)T + B{uj,T)v]} = exp { [^r + (l - e"") v] {(3uj/2y + 0{{(3uj)^)} . 



Moreover (cf. Eqs. §B> and 



A(^) 



26»//32 



l + OiiPco) 



and Eq. ( H6l) yields 



2 fo^sinuz c 
b (2;, t>, r j ~ — / exp I - 



r + 1 - e ^ U; 



(cjV2)}rfcu, 



TT JO Co" 

which coincides exactly with Eq. fl5Bl) . Hence, we have again (see Eqs. f HU]) ^ 



S{z, V, t) ~ Erf 



\{v,t) 



(/?«1), 



(47) 



with A(z, r) defined in Eq. (1371) . In Fig. Hlwe show the goodness of this approximation for 
the small value (3 = 0.1. 

Before proceeding further let us notice an apparent contradiction in the above statements 
on the coincidence of the approximate expressions for the SP between the cases of large 
volatihty and reduced volatility fluctuations, because, at first sight, one would expect the 
opposite, i.e., that high volatility and enhanced volatility fluctuations would be equivalent. 
Let us remind, however, that a large volatility solely means that today's volatility is large 
and this, by any means, does not imply that volatility fluctuations must be intense and vice 
versa. 
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B. Large fluctuations 



We once more start with the exact expression of the SP given in Eq. ( H6i) but supposing 
now intense volatihty fluctuations, so that the vol of vol is greater than the pulling toward 
the normal level. In such 

and from Eq. fl271) we get the approximations 

1 



A(.;) = {[300) 
Hence (cf. Eq. (EH)) 



i/3Lu/2) 



(48) 



(49) 



where, in writing Eq. ( l49l) . we have had to assume that time r is long enough to neglect 
exponential terms of the form e"^'^'^ and higher order terms. This certainly excludes the 
initial stages of the process as will become clearer later. Moreover 

A(lu) 



1 se/zJ^ 









where we have proceeded as in Eq. ( 149|) by neglecting exponentially small terms. Collecting 
these results into Eq. (H6ll we find 



2 f 

S{z,v,t) / 

vr JO 



°° smujz 



exp [- {9t + v) (uj/P)] du, 



(50) 



which after performing the Fourier sine inversion yields 



42| 



71 \9t + V J 

In Fig. IHwe show the accuracy of this approximation for a large value of 
In the original units we have 

2 / S\L — xl \ 

S{x,y,t) arctan — — — — (/3 > 1). 

TT \m'^t + y/a J 



(51) 



(52) 



Let us remark that this expression verifies the boundary condition S{L, y, t) = 0, but not 
the initial condition, since S{x,y,0) 7^ 1 for ?/ 7^ 0. Consequently, Eq. ( l52l) will work only 
after the initial period has elapsed, as otherwise was noted above. We observe, however, 
that the initial condition would be satisfied should we introduce an ad-hoc factor (1 — e"*) 
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FIG. 4: (Color online) Log-log representation of the hitting probability W{z,v,t) = 1 — S{z,v,t) 
as a function of the threshold distance z at time r = 0.5 and v = 9. Diamonds and squares 
represent the numerical evaluation following the exact expression ()30p for two different values of 
f3. The dotted line corresponds to the asymptotic expression for /3 small given by Eq. (j47p . The 
solid line represents the approximate expression obtained trhough Eq. (j5ip . We have chosen the 
realistic parameters of the previous figures and solely modify the parameter k to provide different 
values of /?. 

multiplying the volatility. In this way we obtain a "semi-phenomenological approximation" , 
Sp{x,y,t), which obeys the boundary condition as well as the initial condition: 



where X(t,y) is defined in Eq. ( H2l) . 

V. AVERAGING THE VOLATILITY 

Actual financial data consist in time series of prices from which returns are readily ob- 
tained. Once we have them, volatilities are constructed as the standard deviation of returns. 
Thus, in practice, volatility is a hidden variable that has to be measured by indirect means 




(53) 



19 



fwe refer the reader to Ref. 



44| for more information about this significant question). It 



is, therefore, meaningful to know whether the price of a given asset has reached for the first 
time some critical level regardless the value of its volatility. Let us remark that in physics 



an analogous situation would be knowing whether the position of a Brownian 
first reached some preassigned value without worrying about its velocity 



particle has 
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46|. 



Our goal is then to get the survival and hitting probabilities of the return overlooking 
volatility. We proceed as in Ref. and average the volatility away from S{z,v,t). In 
order to perform such an average we must choose an appropriate probability density for 
it. To this end we make the assumption that the entire process described by Eqs. (EI)-® 
commenced in the infinite past. Consequently, at the present time where we measure the 
return, the volatility has reached the stationary state. We, therefore, define the averaged 
SP, S{z, r), as 

S{z,t)= S{z,v,T)pstiv)dv, (54) 
Jo 

where Pst(f), the stationary probability density of the volatility, is given by the Gamma 
distribution (cf. Eqs. (15]), and (|^): 



Pst[V) 



(55) 



Substituting the exact SP, as given in Eq. (j46l) . into Eq. fl5^ and performing the integral 
over the volatility we have 



S{z,r) 



2 

7r Jo 



_{^Ji+{u) + /i_(^)e-AH-)(l + B{uj, r))_ 

We will write this expression in a more convenient and simplified form. Let us note from 
Eq. ([26]) that 

/i(a;) + [1 + /i_(a;)] + - /i+(u;)]e-^H- 



1 + B{uj,t) 
but 1 + = /i+(co') and 1 — /^+(ci;) 

1 + B{u,t) 

Therefore, 



— /i_(co'). Hence 



/i2(cu)e-^(")^ 



S{z,t) 



2 

TT Jo 



-1 26/ (3^ 



-(2er//32)/._H^^_ 



(56) 
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This is the exact expression for the return's survival probabihty when the volatihty has been 
"thermahzed" . Before proceeding further in seeking approximate expressions for S{z,t), as 
we have done above with S{z,v,t), we will first obtain the SP had the return followed the 
Wiener process. 



A. The Wiener process 

is the geometric Brownian motion which was 



Undoubtedly the most spread market mode 
proposed by Osborne in the late ninety fifties 



471]. In this model the volatility cr is constant 



and the zero-mean return is described by the stochastic differential equation 

dX{t) = adW{t), 

that is, X{t) is the Wiener process with variance o"^. 

Let us denote by So{x,t) the survival probability of the Wiener process inside the semi- 
infinite strip L < X{t) < oo if x > L (or inside — oo < X{t) < L if x < L). This function 
obeys the equation jssl 

dSo_l 2^So 

with initial and boundary conditions 

S'o(x,0) = 1, SoiL,t) = 0. 
Getting along the lines of Sec. Ullwe define the new return variable (cf. Eq. ffTOj) ) 

z = \L — x\ 

for which the hitting problem remains unchanged, except that now the boundary is located 
at z = 0. As before the solution is given in the form of a Fourier integral 

2 r'^ ~ 

So{z,t) = — So{uj,t) sin ujzduj, 
n Jo 

where SQ{u,t) obeys the differential equation 

^ = -(c.aV2)^o, 

with the initial condition 

So{u,0)=V{l/u). 
21 



The solution to this initial- value problem is straightforward and reads 
Hence 



71 Jo 

whence 



S„(M)^Erf(^). (57) 

The comparison of the exact SP for the Wiener return, Eq. fl57|) . with Heston's ap- 
proximate SP given in Eq. (14T!) brings up an interesting detail. Thus for long times, large 
volatility or reduced volatility fluctuations, the approximate SP of the Heston model has the 
same form than that of Wiener SP, changing only the function a^t (which for the Wiener 
case is precisely the return's second moment) by the function m^t + (1 — e~"*)(y/a) for the 
Heston case. Let us observe that this function almost agrees with the return second moment 
which can be shown to be 26 1 



{X\t)) = mh +{y- m") (1 - e-°*)/a. (58) 
B. Asymptotic approximations 

We go back to the exact SP of Heston's return, Eq. (!56|) . Looking at the exponential in 
the right hand side of that equation we easily realize that if the (dimensionless) group of 
parameters given by 

{e/p')T = im/P)H (59) 

is large, we can apply the saddle-point approximation to get an approximate expression of 
the SP. Then proceeding as in Sec. IIII Al and using Eqs. (ETl) and ( IMl) we will have 

^(.,r)^^r^e-(^-/^)-^cia;. 

TT JO UJ 

Whence 

S(..r).Erf(^). (60) 

Figure [5] shows the hitting probability, W{z,t) = 1 — S{z,t), with S{z,t) obtained by 
means of Eq. (I60l) . We compare it with the exact result evaluated through Eq. fl56l) . As 
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FIG. 5: (Color online) Log-log plot of the hitting probability W{z,t) = 1 — S{z,t) as a function 
of z for a fixed time (r = 1, approximately 22 trading days) and for different values of /?. The 
solid line represents the approximation given by Eq. (|60p -which, recall, is independent of /?- while 
dots (i.e., circles, diamonds and squares) represent the exact result. Equation (j60p only works for 
very mild volatility fluctuations (note the clear divergence between the exact probability and its 
approximation for (3 = 1). Parameters are the same than those of Fig. [TJ 

we see in Fig. O the approximation works fairly well for very small values of f5, whereas for 
/? = 1 Eq. ( l60l) significantly deviates from the exact result. We remark that although the 
approximation above is meant for small values of f3 it is itself independent of (3. In Fig. [6] 
we represent the dependence on (3 of the hitting probability and this fact becomes apparent. 
In the original units we have 



Recall that this approximation is valid as long as the group of parameters defined in Eq. 



attains a large value. That is to say, either for long time t, or when the normal level m 
is large, or the (normalized) volatility fluctuation f3 is small. Moreover Eq. (|6T|) exactly 




(61) 
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corresponds to the Wiener SP, Eq. (157|) . by switching the constant volatihty of the latter to 
the normal level. We also note this approximation has the same form than that of Eq. fH7]) 
after setting y = 0. 

Since the approximate expression for S{x,t) given in Eq. fpU]) is the same than that of 
S{z,v = 0, r) given in Eq. fHTl) . the default (or uprising) problem when \L\ oo will be 
given by Eqs. (H3l) and (jUj) when v = 0. That is, 



S{z,r)^l-J- —- (|L|^oo), (62) 



and 



W{z,r)^^j- (1^1 -oo), (63) 



which means an acute exponential decrease of default as \L\ oo. This is to be contrasted 
with a very much slow descent of these probabilities to be discussed next. 



C. High volatility fluctuations 



The other approximation that remains is that of high volatility fluctuations. We start 
from the exact expression of the SP given in Eq. (1561) and if (3 is large we may proceed as 
in Sec. IIV A[ We, therefore, use the expansions given in Eq. ( HHl) along with 



p-liuj) = (J3UJ/2)- 



to write 



/ii(cu) - fil{uj)e-^(^)^ 1 + /3cu/2' 
where we have neglected exponential small terms of the form e~^^'^ . Substituting into Eq. 
f l56|) we obtain 



IT Jo 

Moreover as /5 — > oo we have 

2 1 

l + puj/2_ 

and to the lowest order we get 



2 r°° sinujz 



1 + (3(jj/2 



1 + Puj/2 



1 



S{z,t) 



2 f°° sincjz 



vr Jo LO 



1 + 0(1//?2)1 e-('"/^)"rfcu. 
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FIG. 6: (Color online) Log-log plot of the hitting probability W{z,t) = 1 — S{z,t) as a function 
of (3. Circles and diamonds represent the exact result, obtained through numerical integration of 
Eq. (j56p . for z = 0.01 and two instants of time: r = 0.1 and r = 1. Note that for small values of /3 
the hitting probability is insensitive to /5, in agreement with Eq. (j60p . Moreover, for large values 
of 13 we observe the power-law decay predicted in the asymptotic estimate given by Eq. ([66]). i.e., 
W ~ f3~^. Parameters are the same than those of Fig. [TJ 

Whence 

S{z, r) ~ - arctan f ^ ) , (/3 > 1). (64) 

TT \"'^ J 

The asymptotic SP given by Eq. (1641) has the same functional form than that of Eq. (1511) . 
In fact by neglecting the volatility in Eq. ( 15T1) both approximations coincide. A singular 
characteristic of Eq. flMl) is that it satisfies both boundary and initial conditions, while the 
approximation (ISTj) involving volatility solely obeys the boundary condition. Moreover, and 
contrary to the approximation given by Eq. fl6Tl) . the case of intense volatility fluctuations 
given in Eq. (1641) is not by any means related to the Wiener case discussed in Sec IV A[ 
Therefore, Eq. flM|) . and the same applies to Eq. ([51]), represents a distinct characteristic 
of The Heston model -and, by extension, to any SV model- which has no parallel in the 
Wiener model. 

Figure [7] displays the hitting probability W{z, r) when the SP is given by approximation 
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FIG. 7: (Color online) Log-log plot of the hitting probability W{z, r) = 1 — S{z, r) as a function of z 
with T = 0.5 and for three large values of /?. Solid lines display W{z, r) based on the approximation 
given in Eq. (|64p while dots represent the exact result obtained through the numerical integration 
of Eq. (j56p . Parameters are the same than those of Fig. [TJ 

flM|) along with the exact W obtained through the numerical integration of Eq. fISBl) . The 
curves shown correspond to increasingly large values of jS. Let us note the excellent agree- 
ment between the approximation and the exact result even for moderately high values of 
(3. 

The singular character of Eq. (1641) is best expressed when one considers the possibility 
of default (or uprising) for which \L\ oo. In this case, using the asymptotic estimate 

arctan^ ~ 7r/2 - 1/z + 0{l/z^), 

we get 

9t 

(1^1 -oo), (65) 

and 

(1^1 -oo), (66) 

(see Fig. [6]). Comparing this slow decrease in L with the enhanced exponential fall of the 
previous section (see Eqs. (I62l) - (!63|) ). we see that a boost on volatility fluctuations results 
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FIG. 8: (Color online) Log-log plot of the hitting probability W{z,t) = 1 — S{z,t) as a function 
of z with r = 3 and for large and small volatility fluctuations given, respectively, by /? = 10 and 
f3 = 0.1. Circles and squares represent the exact W{z,t) obtained through Eq. (j56p . Curves 
display W{z,t) based on the approximations given in Eqs. (|60p (dotted line) and (j64p (solid line). 
Parameters are the same than those of Fig. [1] 

in a large increase of risk. 



D. The crossing level and increase of risk 

Figure [8] shows the exact hitting probability for both mild and intense volatility fluctua- 
tions, along with the approximations given by Eqs. fl60l) and fl64|) corresponding, respectively, 



to (3 small and large. Considering 2; = |L| [37| the exact expressions of W intersect at certain 
value Ic of \L\ (in Fig. M, h — 0.336). The intersection marks a turning point of the hit- 
ting problem, for when \L\ < 1^ the probability of reaching \L\ is higher for small volatility 
fluctuations than when these are intense. On the other hand, when \L\ > 1^ the situation 
is reversed, because now the hitting probability corresponding to large values of (3 becomes 
increasingly higher with \L\ than that corresponding to small values of (3. All of this is 
in agreement with the asymptotic estimates given by Eqs. (jHSD and fl^B]) which predict a 
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quadratic exponential decay in |L| for /3 small, and a much slower power-law decay of the 
form 1/|L| when (3 is large. 

We can approximately evaluate the crossing level Ic through the asymptotic forms of the 
SP given in Eqs. flBUl) and flMl) . Indeed, due to the good precision of these approximations 
(see Fig. [H]), a relatively accurate estimation of Ic is obtained from the solution of the 
transcendental equation 

which, in turn, shows that the crossing level depends on volatility fluctuations, /3, and also 
on time and the normal level through the combination 6'r, i.e.. 

In Fig. [9] we represent the numerical solution of Eq. (!67ll in terms of /3, with 6 fixed and 
for three instants of time. We see there that Ic increases with time and also with (3. In fact, 
there seems to be a logarithmic behavior of the crossing level with respect to the volatility 
fluctuations: 

/c~log/3, (68) 

as is shown by the fits in Fig. [91 

The dependence on time of the crossing level is pictured in Fig. [10] where we clearly see 
that Ic increases with time, although now there seems to be a power-law behavior 

Ic ~ {eT)\ (69) 

where 7 depends on (3. 

Let us finally return to the the increase of risk as volatility fluctuations get higher. This 
can be visualized by calculating the ratio 

between the hitting probability W{z,t) (evaluated through the exact SP, Eq. (l56l) ) and 
the hitting probability corresponding to the Wiener process Wo{z,t) = 1 — So{z,t), where 
So{z,t) is given in Eq. (1571) . Note that we use Wo{z,t) to construct the ratio W{z,t) 
because for small volatility fluctuations Heston's and Wiener's SPs share the same form. 
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FIG. 9: (Color online) Numerical solution of Eq. (I67p in terms of (3 with a fixed value of the normal 
level of volatility 9 = 1.92 x 10^'^ (the same of the previous figures) and for r = 0.7, r = 1.3 and 
r = 2.0. Note that for approximately /? > 10, Ic slows down its increase. Solid lines represent a 
logarithmic fit. Parameters are the same than those of Fig. [TJ 

For sufficiently large fluctuations, and according to the above discussion, the ratio W 
should greatly increase as \L\ oo. Indeed, when /3 3> 1 and |L| ^ oo we have from Eqs. 
§^ and dM]) 

W{z,t) ^ {n9r/2y/^e'^'/^'\ (71) 

In Fig. [TT] we represent W{z,t) (evaluated from the exact result) as a function of z (i.e., 
|L|) when /? = 10 and r = 3. As we see there the exponential growth of as |L| increases 
is manifest. 

Figure [TT] also shows that for small values of |L| the ratio is smaller than 1, which implies 
a greater probability of hitting label \L\ when j3 is small than when it is large. However, this 
fact, already observed in Fig. [HI has little consequences on risk control due to the smallness 
of \L\. 

We consider the above findings on the relationship between volatility fiuctuations and risk 
one of the key results of the present work which might have useful practical consequences 
on risk management and control. 
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FIG. 10: (Color online) Numerical solution of Eq. (I67p in terms of Or for /? = 1, /3 = 5 and /? = 10. 
Solid lines are the power law given in Eq. ()69p where the exponent 7 depends on (3 and takes the 
values 7 = 0.3293 ± 0.0051 for /3 = 1, 7 = 0.4212 ± 0.0011 for /3 = 5 and 7 = 0.4358 ± 0.0007 for 
f3 = 10. Parameters are the same than those of Fig. [TJ 

VI. SUMMARY AND CONCLUSIONS 

We have approached the issue of risk evaluation and control as a first-passage problem 
within the assumption of stochastic volatility given by the Heston model. The problem is 
solved when one knows the survival probability, S{x,y,t), to a certain critical level L, of 
the return X{t) starting at X{0) = x (usually x = [37]) with initial volatility Y{0) = y. 
Closely related to the SP it is the hitting probability W{x,y,t) = 1 — S{x,y,t). When the 
critical level, L > 0, is positive W represents the probability of an uprising at time t. On 
the contrary ii L < 0, W gives the probability of having, at time t, a loss quantified by L 
(when L —00 this is the default probability). 

For the Heston model, the SP obeys the backward Fokker-Planck equation, Eq. (jS]), with 
initial and boundary conditions given by Eq. ([9]). The problem have been solved exactly 
by means of the Fourier-sine integral given in Eq. fl29|) . We have then proceeded to get 
handy approximations to the exact SP. Thus, in the asymptotic regimes of either long times 
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FIG. 11: (Color online) Log-log plot of the ratio W{z, r) as a function of z when /3 = 10 and r = 3. 
Note the sharp burst as \L\ increases, in agreement with the approximation given by Eq. (I7ip . 
The inset shows the little decrease of the ratio for small values of |L| (see also Fig. [8|). The turning 
point between these two behaviors is approximately located at \L\ = Ic = 0.336. 

or large initial volatilities, the SP goes to the following Gaussian form which constitutes a 
restatement of the central limit theorem: 

' \L-x 



S{x,y,t) ~ Erf 



(72) 



where X{t,y) is defined in Eq. ( l42l) . This approximation is also valid when the amplitude 
of the volatility fluctuations is small. This amplitude is characterized by the dimension- 
less parameter (3 = k/a formed by the ratio between the vol-of-vol k and the strength of 
the reverting force, toward the normal level, measured by a. A remarkable feature of the 
asymptotic expression, Eq. fl72|) . is that it has the same form had the return followed the 
Wiener process (i.e., constant volatility) instead of the Heston model. 

Other consequences of the asymptotic form above are, for one hand, the non-existence of 
a mean first-passage time, since 

as t oo and the SP does not decays fast enough in order to possess a mean first-passage 
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time. On the other hand, when the volatihty variable y increases the SP decreases as the 
power law: 

S{x,y,t)^y-'/\ 

Obviously these properties are shared by the Wiener process because all of them stem from 
the common form given by Eq. ( 1721) . 

An additional asymptotic form of the survival probability is obtained when one considers 
the fluctuations of the volatility which are characterized by the parameter (3. As mentioned 
above, for mild fluctuations for which /5 ^ 1, the approximate SP is still given by Eq. fl72l) . 
However, when /? ^ 1, i.e., for extreme fluctuations, the asymptotic SP is given by the 
following non- Gaussian form: 

t) ~ - arctan ( ^jf ~ \ (/5 » 1). (73) 

Real financial data consist of time series of prices and the volatility is not directly recorded 
and only observed in an indirect way. This hidden character makes it worth averaging out 
volatility from the expressions of S{x, y, t) and thus solving the hitting problem for the 
return alone. The assumption to be made is that the volatility has reached the stationary 
state characterized by the Gamma distribution. 

Following this way we have obtained the exact expression of S{x,t) given in Eq. (jSSD, 
again in terms of a Fourier integral. As before the averaged SP has two different asymptotic 
forms: a Gaussian approximation for long times or small volatility fluctuations and another 
form, which is non-Gaussian, for intense fluctuations. Thus, when the group of parameters 
{m/f3Yt ^ 1 attains a large value (i.e., for long time t, or large normal level m?, or small 
volatility fluctuations (3) the averaged SP is approximately given by the following Gaussian 
form 

which, curiously enough, has the same form than Eq. (!72|) after setting y = 0. 

The same trend is obtained in the case of large volatility fluctuations, since now the 
following non- Gaussian form emerges 

S{x, t)^- arctan ( ^i^^) (/3 > 1), (75) 
TT \ niH J 

which, again, has the same form of Eq. (175]) after setting y = 0. 
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The most striking difference between these two asymptotic forms of the first-passage 
problem appears when considering the extreme risk of default. In such a case L —oo and 
the hitting probability, W{x,t), corresponding to the Gaussian form (1711) quickly decays by 
following a decreasing quadratic exponential. On the other hand, as volatility fluctuations 
increase, the non-Gaussian approximation, Eq. (175 p . decays much slowly by following a 
power law. This means a significant increase of risk as the fluctuations of the volatility soar. 
We have shown that there is a crossing level Ic = hiP, m^t) from which risk exponentially 
rises. For a fixed value of time, the crossing level seems to follow the logarithmic law 
Ic ~ log /3, while for a fixed value of (3 the crossing level apparently increases as a power law 
Ic ~ {mHy, where the exponent 7 depends of p. 

Practical consequences of these findings on actual markets will be presented in forthcom- 
ing works. 
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APPENDIX A: FUNCTIONS A{u;,t) AND B{uj,t) 
Function B{uj,t) is the solution of the Riccati equation 



B = -B-B'^ + {(3u/2f, 



(Al) 



with initial condition B{uj,Q) = 0. Define a new function Z{u!,t) related to B by 



B 



Z 



then Z obeys the linear equation 



Z + Z- {Puj/2fZ = 0, 



whose solution is 



Z{uj,t) = Ci(cj)e'^-(^)" + C2{uj)e 
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where 6*1(0;) and 02(00) are arbitrary and /U±(ci;) are defined in Eq. (1271) . The expression of 
B{ijJ,t) is thus given by 

^ ' ^ l + [C2(a;)/Ci(a;)]e-AH- 

where A(ct;) is defined in Eq. (1271) . From the initial condition we get 



Therefore, 



1 _ Q-M'^h 



(A2) 



l + [/i_(^)//i+M]e-^(^)^' 
which proves Eq. fl26l) . 

We now substitute Eq. (jM]) into Eq. ([231) and define C = e-^^"") as a new integration 

variable, we get 



1-C 



/3V A(^) 7e-AM^ C[l + (/i-^+)C] 



Taking into account 



1-C 



:dC = InC - (1 + /i+//i_) ln[l + (/i+//i-)C], 



C[i + (/i_//i+)C] 

and recalling that H-iuj) + n+iyj) = A{u) (cf. Eq. (1271) ). we finally obtain 



/i+(cu) + /i_(w)e-^(")" 



AM 



(A3) 



which proves Eq. 
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